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METHOD FOR ESTIMATING MINIMUM REQUIRED EJECTION 
VELOCITY FOR PARACHUTE DEPLOYMENT 

By Earle K. Huckins III and Lament R. Poole 
Langley Research Center 

SUMMARY 

A technique was developed for estimating, for preliminary design purposes, the min- 
imum required ejection velocity for forced ejection parachute deployments. The approach 
used was to obtain a first integral to the equation for the linear motion of an inelastic para- 
chute deployed in the lines-first mode. A criterion was developed which limits the theory 
to forced ejection deployments by constraining the ratio of acceleration of the deployment 
bag to acceleration of the vehicle at bag strip. In order to perform the integration, the 
acceleration of the deployment bag relative to the vehicle was approximated by a function 
of deployed distance. The theory was extended for application to static vertical ground 
tests and cases in which the bag-strip velocity is constrained to be a specified percentage 
of the ejection velocity. The results obtained by using the present approximation closely 
agreed with results obtained by numerical integration for a series of hypothetical sample 
cases. 


INTRODUCTION 

A method of lines-first (ref. 1) parachute deployment, which is particularly useful 
for parachute deployment in the wake of a bluff body having a low ballistic coefficient, is 
forceful ejection of the parachute pack from the vehicle. A fundamental design parameter 
for the associated deployment system is the pack ejection velocity. In the past, selection 
of the design ejection velocity has been based primarily on previous experience with para- 
chutes of similar size and weight. 

If, as in many cases, the time interval of deployment is not constrained, minimizing 
the design ejection velocity is usually desirable. A minimum ejection velocity results in 
minimum reaction force acting on the vehicle, since the reaction load for an ideal deploy- 
ment system with a fixed stroke is directly proportional to the square of the ejection 
velocity (ref. 2). In addition, a reduction in the reaction load is usually accompanied by 
a reduction in system weight, both of which are desirable system characteristics. 

A definition of the minimum required ejection velocity must be based on some arbi- 
trary criterion for successful deployment. The simplest definition, and the one which is 



used in the present analysis, is that the minimum required ejection velocity is the lowest 
ejection velocity for which the parachute canopy will be completely extracted ("stripped") 
from the deployment bag. Under this criterion, the deployment will be satisfactory in the 
sense that the parachute will be completely extended and exposed to the flow field. How- 
ever, satisfactory inflation and proper functioning of the parachute depend on numerous 
additional factors. No experimental evidence is presently available to verify or refute 
the assumption that complete bag strip corresponds, in itself, to successful deployment. 
However, defining the required ejection velocity as that required to strip the deployment 
bag completely, produces some results which are quite useful in preliminary system 
design. 

The appropriate pack ejection velocity for any particular set of constraints can be 
accurately determined by numerical integration of the equations for the motion of the 
deployment bag relative to the vehicle (ref. 2). However, numerical integration techniques 
are not always practical during the preliminary design phase. In the present analysis, a 
closed-form approximation to the minimum required ejection velocity (as defined above) 
is developed which accounts for the dissipative influences of vehicle drag and bag friction, 
in addition to the effect of deployment-bag drag in an arbitrary wake profile. Criteria 
are developed which bound the range of validity for the approximate result. 

SYMBOLS 


A reference area, meters^ 

Drag 

Cj) drag coefficient, 

q A 

oo 

Dq nominal diameter of parachute, meters 

Fj,g unfurling resistance force, newtons 

g acceleration due to gravity, meters/second^ 

K integral parameters defined by equations (11) to (13), meters/kilogram 

I extended length of segment of parachute, meters 

m mass, kilograms 

m' linear mass density of parachute exiting deployment- bag mouth, 

kilograms/meter 


2 



free-stream dynamic pressure, newtons/meter 2 
Tg tension at mouth of deployment bag, newtons 

t time, seconds 

V velocity vector, meters/second 

X deployed distance, meters (see fig. 1) 

^ ballistic coefficient, kilograms/meter2 

y vehicle flight-path angle, degrees 

e ratio of bag-strip velocity to ejection velocity 

ri wake parameter (ratio of drag coefficient of deployment bag in vehicle wake 

to that in free stream) 

(AV)^gq required ejection velocity corresponding to a specified velocity at bag strip, 
meters/second 

Subscripts: 

b deployment bag 

bs bag strip 

c parachute canopy 

cs canopy skirt 

e deployment bag plus contents 

p parachute 

si suspension lines 

T total parachute and deployment bag 

3 


I 


u 


unfurled portion of parachute 


V vehicle 

w vehicle wake 

1 interval of deployed distance from x = 0 to x = Zgi, meters 

2 interval of deployed distance from x = Zgi to x = Zg^ + meters 

Dots and bars over symbols denote time derivatives and average values, respectively. 

ANALYSIS 

In the analysis to follow, an approximate first integral (relative energy integral) to 
the equation which governs the acceleration of the deployment bag relative to the vehicle 
is obtained. This integral can be used to determine the ejection velocity which corre- 
sponds to a specified relative velocity at bag strip. 

Equation of Motion 

A typical vehicle-parachute configuration for a deployment of the lines-first type is 
shown in figure 1 . The forces which influence the dynamics of this type of deployment 
process are shown in figure 2. An equation for the acceleration of the deployment bag 
relative to the vehicle was derived in reference 2 under the following assumptions: 

(1) The motion is linear. 

(2) The parachute is inelastic. 

(3) The partially unfurled parachute is in tension during the deployment process. 

(4) The aerodynamic drag of the partially unfurled parachute is much less than the 
drag of the vehicle. 

(5) The deployment rate is much less than the velocity of the vehicle. This assump- 
tion allows the dependence of the dynamic pressure at the deployment bag on the deploy- 
ment rate to be neglected. 

Under these assumptions, the governing equation of relative motion was derived in 
reference 2 and can be rewritten in the following form: 

- I're 

X 

mg m.^^ + my 


4 





where 


X deployment acceleration 

ri wake parameter 

Cjj drag coefficient 

A reference area 

free-stream dynamic pressure 
Fre unfurling resistance force 

m mass 

Tg tension at mouth of deployment bag, m'x^ + Fj.g 

m' linear mass density of parachute exiting deployment-b^ mouth 

X deployment velocity 

and subscripts: 

b deployment bag 

e deployment bag plus contents 

V vehicle 

u unfurled portion of parachute 

In deriving the governing equation, the effect of premature canopy inflation on the 
motion of the deployment bag has been neglected. Limited canopy inflation prior to bag 
strip results in a wedging effect on the bag mouth which tends to assist stripping of the 
deployment bag. However, this process is quite difficult to model mathematically and 
strongly depends on factors, such as packing irregularities, that are random in nature. 
The effect of neglecting this influence is to generate results which are conservative by 
a small but indeterminate amount. 
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Approximate First Integral 


The velocity of the deployment bag relative to the vehicle can be determined as a 
function of deployed distance by using the procedure given below. The relative accelera- 
tion can be expanded by the chain rule for differentiation as 


r. _ dx _ dx dx _ dx 
dt ~ dx dt “ dx 


( 2 ) 


Therefore, 

X dx = X dx 


(3) 


which can be integrated to give 



sl'^^c 




where 


(4) 


Zgj extended length of suspension lines 

extended length of canopy 

x^jg relative velocity at bag strip (i.e., at x = Zg^ + l^) 

The above integration assumes implicitly that the deployed distance is a monotonically 
increasing function of time. For a zero relative velocity at bag strip, a necessary condi- 
tion for monotonic increase in the deployed distance is a nonpositive acceleration of the 
deployment bag at bag strip. By use of equation (1), this condition can be expressed as 
follows : 

^(CpA)^q^ - 

”^b < 1 (5) 

(CpA)^q^ + Fre 

+ mp 

where mp is the total mass of the parachute and all parameters are evaluated at the 
time of bag strip. This criterion limits the theory to forced ejection deployments. For 
cases in which the criterion is violated, the deployment process approaches that of a 
drogue deployment and the bag- strip velocity is always greater than zero. 

The required ejection velocity is equal to the deployment velocity evaluated at the 
lower limit (x = 0). Solving equation (4) for the required ejection velocity gives 
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( 6 ) 


where 



I 


«=^sl+^c .. ^ 

X dx 

x=0 


pX=i 

Jx=0 


-Zgl+^c 


7?(CDA)^^q^ - Fre (CpA)^q 


+ Ttj 

V OO -D 


nif 


m-,r + m, 


dx 


(7) 


In order to obtain the integral of the relative acceleration, the relative acceleration 
must be approximated by a function of the deployed distance only. This approximation 
should be conservative with regard to ultimate determination of the required ejection 
velocity. A convenient approximate expression for the relative acceleration can be 
obtained by introducing the following assumptions: 

(1) The wake parameter and the unfurling resistance are functions of the deployed 
distance. In addition, the unfurling resistance has a constant value during deployment of 
the suspension lines and another constant value during deployment of the canopy. 

(2) Drag coefficients are constant. 

(3) The mass of the unfurled portion of the parachute is negligible in comparison 
with the mass of the vehicle. 

(4) The tension at the deployment-bag mouth can be conservatively approximated by 
assuming a constant deployment rate equal to the ejection velocity and by neglecting the 
unfurling resistance. These assumptions are used only in approximating the tension. For 
a case in which the predicted required ejection velocity is relatively high, this assumption 
makes the result highly conservative. 

(5) The free-stream dynamic pressure is constant during the deployment process. 
This assumption is conservative for a negative dynamic pressure gradient. For a con- 
stant free-stream density during the deployment process, a negative dynamic pressure 
gradient results under the following condition: 


^^g_siiM^_^ 1 
m’x2 + Fre ^ ~ 

where jSy is the ballistic coefficient of the vehicle. Integrations for cases 
this criterion are valid but the results must be considered nonconservative. 


(8) 


which violate 


Under the above assumptions, the integral of the relative acceleration (eq. (7)) is 
given by 
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( 9 ) 


j dx = (coAU -L dx 

Jx=0 ' /b °oJv=n “p 


S 

x=0 “"e 


(^-)l C'‘ ^ *■ ■ Cf ^ ^ 

ol 

^?sl + k) - s5(%=o)^ 


where 


(Fre)j unfurling resistance over interval of deployed distance from x = 0 to 

X = Isi 

(Fre )2 unfurling resistance over interval of deployed distance from x = Zgj to 
^ ” ^sl ^c 


An expression for estimating the required ejection velocity can now be generated 
by substituting equation (9) into equation (6). The final result can be written as 

fq 

f{lBl * Ic) * Ki(Fre)i ^ K2(Fre)2 - / 


(AV) 


1 \V2 C 


req 


,1 - 2 




where 



mv/ V. 

-j 

(10) 

Ki = 

‘^x=0 ® 

(11) 

K2 = 


(12) 

Kw = 

r^^sl+^c jLdx 

‘^x=0 *^e 

(13) 


The integral parameters defined by equations (11) to (13) can be evaluated numerically, 
analytically, or graphically for a particular parachute and vehicle wake profile. 

For some applications, it may be useful to determine the ejection velocity which 
corresponds to a bag-strip velocity that is a specified fraction of the ejection velocity. 
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That is, 


^bs ~ 


(14) 


Under this criterion, equation (10) reduces to 


(AV) 


req 


1/2 


(AV) 


^bs“^(^V)j.0q 


req 


1-2 


%s" 


(15) 


For small values of e and the ratio of parachute mass to vehicle mass, equation (15) 
can be closely approximated by 


(AV) 


req 


^bs=^(^V) 


- (l + 2 ^^)(^V)reql 


req 


^bs=0 


(16) 


The significance of this equation is to indicate the high sensitivity of bag-strip velocity to 
ejection velocity for ejection velocities close to the minimum required. 

An expression for the required ejection velocity for a static vertical ground test, 
which corresponds to equation (10) for an in-flight deployment, can be derived in a simi- 
lar manner using equation (1). In a ground environment the effective deceleration of the 
vehicle, due to drag and tension (the second term of eq. (1)), is the local acceleration due 
to gravity. The deployment-bag drag acts in a direction opposite to the relative deploy- 
ment motion. The wake integral parameter is equal to the sum of the other two integral 
parameters. The free-stream dynamic pressure varies considerably between pack ejec- 
tion and bag strip. Therefore, an average dynamic pressure must be estimated in order 
to perform the closed-form integration. For these conditions, the required ejection 
velocity is given by 


(AV) 


req < bs 


= + 2 




S(^sl ^c) + ^l(^re)i + ^2(^re)2 

1/2 


(17) 


where q^ is the average dynamic pressure over the appropriate interval of the unfurling 
process. 


Evaluation of Integral Parameters 


The integral parameters defined by equations (11) to (13) express essentially the 
sensitivity of the required ejection velocity to the mass distribution of the parachute. 
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These integrals can be evaluated graphically or numerically for a particular parachute 
mass distribution and an experimentally determined wake profile. In many cases, the 
parachute mass distribution can be closely approximated by that of an equivalent flat cir- 
cular parachute. Closed-form integrations for an ideal flat circular parachute are devel- 
oped in the appendix and can be expressed as 

(18) 


(19) 

where 


1 nisi ^ r mrrJ 


and 


K 2 = - tan- 

2^/mcmb \"^b/ 


1/2 


m mass 

Dq nominal diameter of parachute canopy 

and subscripts: 
si suspension line 

T total parachute and deployment bag 

b deployment bag 

c parachute canopy 

The form of the integral parameter K 2 given in equation (19) can be used to study 
the effects of variations in bag weight. Increasing the bag weight by ballasting minimizes 
the relative deceleration due to unfurling resistance and decreases the relative accelera- 
tion due to bag drag. However, ballasting may in some cases be undesirable because of 
the associated weight penalty. 

The integral defined by equation (13) can be approximated by 
^w ~ ^1^1 

where rj is an average value of the wake parameter over the appropriate interval. 


I 
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RESULTS AND DISCUSSION 


In order to study the accuracy obtained by using the present approximation to the 
required ejection velocity, a series of hypothetical parachute-deployment situations were 
studied. These sample cases were chosen such that the criteria given by equation (5) and 
equation (8) were satisfied. An ideal flat circular parachute was assumed. The effects of 
variations in parachute size and weight, ballistic coefficient of the vehicle, and free-stream 
dynamic pressure on prediction accuracy were studied. Parameters for the selected test 
cases are given in table I. The required ejection velocity was calculated from equa- 
tion (10); it was determined for the sample cases by numerical integration of equation (1) 
with a digital computer. Results obtained by using both techniques are given in table II. 
Agreement between the two methods of calculation is good. The largest conservative 
error (11%) is found in case 3. This particular example postulated a relatively large 
parachute deployed from a vehicle having a low ballistic coefficient at a high dynamic 
pressure. For these conditions, a relatively large required ejection velocity would be 
expected. The assumption of tension based on a constant deployment rate results in an 
overprediction when equation (10) is used. However, the result is not considered to be 
overly conservative and could be used for engineering design purposes. 

CONCLUDING REMARKS 

An analysis has been made to determine the minimum required ejection velocity for 
parachute deployment by forced ejection. An approximate first integral to the equations 
for the inelastic linear motion characteristic of a deployment of the lines-first type was 
obtained. Criteria were developed which insure conservative results and limit the theory 
to forced ejection deployments. Results obtained by using the approximate closed-form 
integration closely agreed with results obtained by numerical integration. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 30, 1971. 


13 


I 


APPENDIX 


CALCULATION OF APPROXIMATE INTEGRALS 

Approximate values of the integrals describing the effects of unfurling resistance 
and the wake of the vehicle on the required ejection velocity are derived in this appendix. 
In the derivation, it is assumed that the parachute consists of an ideal flat circular canopy 
and uniform suspension lines. A mass distribution for this type of parachute is shown 
in the following sketch: 



X 


Sketch 1 

where 

m’ linear mass density 

I unfurled length 

and subscripts: 
si suspension lines 
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APPENDIX - Continued 


cs 


canopy 


canopy skirt 


The unfurling resistance parameter for the suspension lines can be written as 
follows: 




(Al) 


where 


me(x) 


displacement of deployment bag relative to towing vehicle 

mass of deployment bag and its instantaneous contents as a function of 
deployed distance 


In general, mg(x) can be represented by 
mg(x) = mrp - y m'dx 


(A2) 


where mrp is the mass of the entire parachute and the deployment bag. In the range of 
the suspension lines, equation (A2) can be written as 

mg(x) = mrp - mgjX (A3) 

By substituting equation (A3) into equation (Al) and integrating, the parameter Kj 
can be written as 


Ki = - ln(l - 
1 mgi \ m-r; 


(A4) 


where mgj^ is the mass of the suspension lines. 

The unfurling resistance parameter for the canopy can be written in general form as 


„ _ r'si+'c 1 

By defining Xg = x - Zgj^ the parameter can be written as 
’'h 


(A5) 


K 2 =r-^— dx, 

Jo mg(Xg) 


(A6) 
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APPENDIX - Concluded 


In general, can be written as 


r^c 

me(Xc) = dx^ 

By referring to sketch 1, nic(Xg) can be written as 


m 




By substituting equation (A8) into equation (A7) and integrating, 


(A7) 


(A8) 


^e{^c) = niT - mgi 


"^cs^'c 



(A9) 


where Dq is the nominal diameter and is equal to Through substitution of equa- 

tion (A9) into equation (A6) and integration, there results 



(AlO) 


where 


m^ mass of parachute canopy 

m ]3 mass of deployment bag 


The integral parameter which describes vehicle wake sensitivity can be written in 
the following form : 


Kw 


r^si'^^c 

^0 


V 

me(x) 


dx 


(All) 


where t? is a ratio of the deployment-bag drag in the vehicle wake to the deployment- bag 
drag in the free stream. 

The wake integral parameter can be approximated by 

Kw « rjjKi + lygKg (A12) 

where rj is the average value of 77 over the appropriate interval of deployed distance. 
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TABLE I.- PARAMETERS USED IN SAMPLE CALCULATIONS 


Vehicle characteristics: 

Weight = 5000 newtons 
Drag coefficient = 0.5 
Wake parameters: 

Vi = 0.4; 772 

Parachute characteristics: 

Ideal flat circular canopy 

^sl " 

Ic = Do/2 

msi = (m^ - mjjj/3 
(Fre)j = 50 newtons 

(Fre )2 = 200 newtons 
Bag weight = 10 newtons 

^Cj)A^^ = 0.1 meter 2 

Conditions at parachute pack ejection: 
Altitude = 20 kilometers 
Flight path range = -90° 



^sl ^c’ 
m 


15 

15 

30 

30 


m 

kg/m 2 

50 

200 

50 

200 


^CO’ 

N/m2 

3000 

2000 

3000 

2000 
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TABLE n.- RESULTS OF SAMPLE CALCULATIONS 


Case 

Integral parameters, 
m/kg 

m/sec 

Kl 

K2 

Kw 

Approximate 
method 
(eq. (10)) 

Numerical 
integration 
(eq. (D) 

1 

0.5619 

1.7500 

1.6248 

41.46 

40.33 

2 

.5619 

1.7500 

1.6248 

21.01 

21.64 

3 

.4656 

2.3764 

2.0873 

64.66 

58.16 

4 

.4656 

2.3764 

2.0873 

30.86 

30.27 


NASA-Langley, 1971 31 


L-7704 


19 



